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We introduce a new cosmological diagnostic pair {r, s} called Statefinder. The Statefinder is 
a geometrical diagnostic and allows us to characterize the properties of dark energy in a model 
independent manner. The Statefinder is dimensionless and is constructed from the scale factor of 
the Universe and its time derivatives only. The parameter r forms the next step in the hierarchy of 
geometrical cosmological parameters after the Hubble parameter H and the deceleration parameter 
q, while s is a linear combination of q and r chosen in such a way that it does not depend upon the 
dark energy density. The Statefinder pair {r, s} is algebraically related to the equation of state of 
dark energy and its first time derivative. The Statefinder pair is calculated for a number of existing 
models of dark energy having both constant and variable w. For the case of a cosmological constant 
the Statefinder acquires a particularly simple form. We demonstrate that the Statefinder diagnostic 
can effectively differentiate between different forms of dark energy. We also show that the mean 
„Cj ■ Statefinder pair can be determined to very high accuracy from a SNAP- type experiment. 
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Recent observations of type la supernovae indicate that the expansion of the Universe is accelerating rather than 
, slowing down [1]. These results, when combined with cosmic microwave background (CMB) observations of a peak in 
J> ■ the angular power spectrum on degree scales [2,3], strongly suggest that the Universe is spatially flat with ^1/3 of the 
OO \ critical energy density being in non-relativistic matter and ~ 2/3 in a smooth component with large negative pressure 
('dark energy'). Indirect support for dark energy (known long ago) comes from the examination of gravitational 
clustering within the framework of the standard gravitational instability scenario (see the reviews [4,5]). Finally, 
with recent data on the galaxy power spectrum from 2dF Galaxy Survey combined with CMB data, the existence of 
' dark energy can be proved without using the supernovae data at all [6]. A large body of recent work has focussed 
on understanding the nature of dark energy and its possible relation to a fundamental theory of matter such as 
***^; M-theory, supergravity etc. Despite the considerable effort in this direction, both the nature of dark energy as well 
'■^ | as its cosmological origin remain enigmatic at present. 

The simplest model for dark energy is a cosmological constant A, whose energy density remains constant with time 
e = A/8irG and whose effective equation of state remains fixed, w = P/e = — 1 (P is the pressure) as the Universe 
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evolves. The cold dark matter (CDM) model with the cosmological constant having the corresponding mass density 
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where h is the Hubble constant Hq in terms of 100 km s _1 Mpc _1 and ft a = 0.7 ± 0.1, h — 0.7 ± 0.1, provides an 
excellent explanation for the acceleration of the Universe and other existing observational data. However, it remains 
quite possible that the dark energy density may depend sufficiently weakly upon time. This follows from many 
proposed models. The possibility that dark energy could be dynamical is also suggested by the remarkable qualitative 
analogy between the observed properties of dark energy and properties of a different type of 'dark energy' - namely 
the inflaton field - postulated in the inflationary scenario of the early Universe. 

Once we allow the dark energy density to be time-dependent, then the next simplest class of models are those with 
a constant, non-positive w. We shall call this class 'quiessence' (Q) (w < —1/3 is a necessary condition to make the 
universe accelerate). Examples include a tangled and 'frustrated' network of cosmic strings w = —1/3 and domain 
walls w = —2/3. More generally, in a Friedmann-Robertson- Walker (FRW) background with the presence of CDM, 
an arbitrary but constant w for dark energy from the range (—1,0) can be achieved by using a scalar field with a 
hyperbolic sine potential (see Eq. (9) below) [4] . It may be noted that in principle the value of w may be even less 
than —1; the present observational data do not exclude this possibility but limit the constant w in the range of about 
(-1.6,-0.8) [7]. 

A more generic alternative to A and Q is presented by 'kinessence' (K) which refers to dark energy with a time 
dependent w. Examples of kinessence include 'quintessence' - a scalar field <p with a self-interaction potential V(<j>) 
minimally coupled to gravity (see [4] for numerous references), as well as the 'Chaplygin gas' model [8] and braneworld 
models of dark energy [9,10]. These three alternatives are summarized in the Table I (where z = a(to)/a(t) — 1 is the 
redshift, a(t) is a FRW scale factor and the subscript denotes the present moment). 
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Dark Energy 



State Parameter 



Energy Density Parameter 



Cosmological constant 
Quiessence 

Kinessence 



w(z) = constant = — 1 
w{z) = constant < — 1/3 
w(z) / constant 



p{z) — A/8tyG — constant 
e(z) =e (l + z) 3(1+w) 



e(z) — eo exp 



3/«> 



1+z' 



TABLE I. 



The effective equation of state is clearly an important property of dark energy. This has led to numerous attempts to 
reconstruct the former from observations of high redshift supernovae in a model independent manner [11-13]. However, 
for field-theoretical models of dark energy, the equation of state is not a fundamental property. Strictly speaking it has 
reference only to an exactly isotropic FRW background. For small perturbations superimposed on a FRW background, 
the pressure tensor is generically non-diagonal (non-barotropic) , and the velocity of signal propagation need not be 
given by the standard hydrodynamic expression yJdP/de. Moreover, the very notions of e and P for dark energy 
pre-suppose the Einstein interpretation of gravitational filed equations (not to be confused with the notion of the 
Einstein frame which is used in scalar-tensor and string theories of gravity!). Namely, even if the real equations for a 
given model are not the 4-D Einstein equations at all (examples include dark energy models in scalar-tensor [14] and 
brane [9,10] gravity), one can still write them formally in the Einstein form, by placing the Einstein tensor Rij — 
into the left-hand side, and by grouping all other terms in the right-hand side and calling them (after dividing by 
SttG) 'the effective energy-momentum tensor of matter'. After that, the energy-momentum tensor of dust-like matter 
(describing CDM and baryons) is subtracted from the latter, and the remaining part is used to define e and P for 
'dark energy'. All this reveals how ambiguous the notion of 'equation of state' can be for a non-Einsteinian model of 
dark energy. 

Fundamental variables (at least, at the field-theoretical level of consideration) are either geometrical (astronomical) 
if they are constructed from a space-time metric directly, or physical - those which depend upon properties of 
physical fields carrying dark energy. Physical variables are, of course, model-dependent, while geometrical variables 
are more universal. Additionally, the latter do not depend upon uncertainly measured physical quantities such as the 
present density of dust-like matter fl m . That is why we emphasise the use of geometrical variables when describing 
the present expansion of the Universe and properties of 'dark energy'. 

The oldest and most well-known geometric variables are the Hubble constant Ho and the current value of the 
deceleration parameter qo. At present, accurate measurements of the expansion law of the Universe during the 
past are also possible (e.g., using the luminosity distance to distant supernovae), therefore these variables should 
be generalized to the Hubble parameter H(t) = a/a and the deceleration parameter q(t) = —ad/ a 2 = —d/aH 2 
(H = H(t ) and q = q(to)). However, both the necessity of consideration of more general models of dark energy 
than a cosmological constant, and the remarkable increase in the accuracy of cosmological observational data during 
the last few years, compel us to advance beyond these two important quantities. For this reason, in this letter we 
propose a new geometrical diagnostic pair for dark energy. This diagnostic is constructed from the a(t) and its 
derivatives up to the third order. Namely, we introduce the Statefinder pair {r, s}: 
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r(z) is a natural next step beyond H(z) and q(z). We will soon see that it has a remarkable property for the basic flat 
ACDM FRW cosmological model. s(z) is a linear combination of r(z) and q(z). In a companion paper we shall show 
that a particular combination of two variables from the above three e.g., q and s, can provide an excellent diagnostic 
for describing the properties of dark energy [15]. 

Below we will assume that the Universe is spatially flat, k = 0. This assumption naturally follows from the 
simplest versions of the inflationary scenario and is convincingly confirmed by recent CMB experiments [3]. At 
late times (z < 10 4 ) the Universe is well described by a two component fluid consisting of non-relativistic matter 
(CDM+baryons) Q m and dark energy Ux = 1 — TO . In this case the Statefinder pair acquires the form 
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r = l + -n x w{l + w)--n x -, (3) 

S = 1 + W ~3^H (4) 

where w = Px/sx- Thus, if the role of dark energy is played by a cosmological constant (w = —1), then the value of 
r stays pegged at r = 1 throughout the entire matter dominated epoch and at all future times; i.e., r = 1 for z < 10 4 
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irrespective of the current value of f2 TO . The extreme simplicity of the parameter r(z) for the basic cosmological 
model (ACDM) which also provides the best fit to existing observational data may, in fact, prove not to be a mere 
coincidence ! 1 Very different behaviour is predicted for quiessence and kinessence for which r is a function of time. 
In particular, if dark energy is attributed to a minimally coupled scalar field <j> (quintessence), 

r=1 + i^ + ^. <„ 

The properties of the second Statefinder V complement those of the first. For the basic ACDM model with any 
non-zero A, s = 0. Moreover, s does not depend neither on time, nor on Q m , for quiessence models for which s = 1 + w. 
In marked contrast, s generically depends on time for kinessence. E.g., for quintessence: 

cf) 2 -2V 

Thus, properties of the Statefinder pair {r, s} enable it to differentiate between the three canonical forms of dark 
energy described in Table 1. 

It is straightforward to invert Eqs. (3,4) and express w and w in terms of the Statefinder pair. However, w is more 
directly related to the deceleration parameter: 

Thus, w a composite quantity since it is constructed out of physical (f2jc) as well as geometrical (q) parameters. Note 
that for quintessence, w > — 1 but w may have any sign (for models with w > and w < the epoch of dark energy 
domination is usually a transient). The relationship between geometrical and physical parameters is summarized in 
Table 2. 



TABLE II. Relationship between geometrical and physical parameters characterizing the observable Universe 



Geometrical parameters Related physical parameters 



H = a/a 


^total , ^curvature 


q — —a/aH 2 
r = a/aH 3 


tli,Wi 


£li,Wi,ibi 


s = (r - l)/3(q - 1/2) 


Wi,Wi 



Let us now study the Statefinder pair for different models of dark energy in greater detail. As was mentioned 
already, its value is equal to {1,0} for any ACDM model with a non-zero A. Quiessence models (QCDM) have a 
constant w, asa result 
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r = 1 + -fl Q w(l + w), s = l + w. (8) 

Two values of the equation of state are singled out for special attention: w = —1/3 (cosmic strings) and w = —2/3 
(domain walls). In both cases the first Statefinder has the simple form r(t) = 1 — fig(f) = fi m (t). As a result, r(t) — > 1 
for t <C to, r(t) — > for t 3> t and r ~ 0.3 at the present time when fig (to) — 0.7. This leads to a degeneracy in 
ro for the dual value w = —1/3,-2/3. Though generic, this degeneracy is easily broken when one adds information 
from the second Statefinder s. Note that the case of an arbitrary — 1 < w < in the presence of a non-zero fi m can 
be achieved using quintessence with the potential [4] (see also [17]) 2 



1 Note that the quantity r(z) was also considered in the paper [16] for a non-flat case when it is time-dependent. However, its 
remarkable property for the flat ACDM model was not emphasized. For completeness, let us mention that r — 2q = Q m (z), s = 
2/3 for a matter dominated non-flat CDM model with negligible amounts of dark energy and radiation. 

2 There are some misprints in numerical coefficients in Eqs (119-121) of [4] which are corrected here. 
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In this case, r < 1, < s < 1. 

Let us now turn to the quintessence case where r and s are given by Eqs. (5) and (6) correspondingly. To this 
category belong scalar fields with 'tracker' potentials, for which the scalar field <j> approaches a common evolutionary 
path from a wide range of initial conditions [18]. Tracker potentials satisfy V"V/(V') 2 > 1. We consider the simplest 
case of an inverse power-law potential V{<p) — V~o/(f> a , a > first studied in [19]. For this potential, the region of 
initial conditions for <j> for which the tracker regime has been reached before the end of the matter-dominated stage 
is (f>i n <C Mp = \/G, and the present value of quintessence is 4>(to) ~ Mp. The evolving values of the Statefinder 
pair for this potential with a = 2 and a — 4 are shown in Figs. 1 and 2. Also shown are results for the cosmological 
constant and quiessence. During tracking ^- oc t 4 /( 2+Q ) as a result quintessence always becomes dominant at late 
times. The equation of state of quintessence and the corresponding value of the Statefinder pair is given by 

w= — , r w 1, s w 1 + w (10) 

a + 2 

(wb = 1/3,0 during the radiation- and matter-dominated epochs respectively). 
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FIG. 1. The Statefinder pair (r, s) is shown for different forms of dark energy. In quiessence (Q) models (w =constant ^ — 1) 
the value of s remains fixed at s = 1 + w while the value of r asymptotically declines to r(t 2> to) ~ 1 + + w). Two 
models of quiessence corresponding to wq — —0.25, —0.5 are shown. Kinessence (K) models are presented by a scalar field 
(quintessence) rolling down the potential V(cj>) oc 4>~ a with a = 2,4. These models commence their evolution on a tracker 
trajectory described by (10) and asymptotically approach ACDM at late times. ACDM (r — l,s = 0) and SCDM in the 
absence of A (r = 1, s = 1) are the fixed points of the system. The hatched region is disallowed in quiessence models and in 
the kinessence model which we consider. The filled circles show the current values of the Statefinder pair (r, s) for the Q and 
K models (fi m o = 0.3). 



Constraints from structure formation and the CMB suggest that dark energy must be subdominant at z > 1. 
Primordial nucleosynthesis arguments impose the stringent constraint: fix < 0.05 at z ~ 10 9 [20]. Small values of 
fix and w substantially decrease the terms flxw and Qxw/H which appear in the RHS of (3) and ensure that the 
Statefinder r remains close to unity at high z. This is exactly what one finds from Fig. 2. The extreme sensitivity of 
r to an evolving equation of state of the tracker field is reflected by the fact that the value of r declines rapidly as the 
Universe expands, dropping to ~ 50% from its starting value by z ~ 1, even though dark energy remains subdominant 
at this epoch. 

As is apparent from Fig. 2, the discriminating power of r and s can be significant even at moderate redshifts. 
Since and CIq usually decrease faster with redshift than fix, the value of r(z) for both the cosmological constant 
and quiessence is generally closer to unity at a given large redshift than the corresponding value for a tracker field 
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FIG. 2. The Statefinder pair {r, s} is shown for dark energy consisting of a cosmological constant A, quiessence 'Q' with an 
unevolving equation of state w — —0.8 and the inverse power law tracker model V = Vo/<j> 2 , referred here as kinessence 'K'. 
The lower left panel shows r(z) while the lower right panel shows s(z). Kinessence has a time-dependent equation of state 
which is shown in the top right panel. The fractional density in matter and kinessence is shown in the top left panel. 



(kinessence). Thus, whereas the current value of r$ allows us to differentiate A from Q and K, the value of r at 
moderate redshifts distinguishes K from A and Q. This feature is even more pronounced in the second Statefinder s, 
whose value does not explicitly depend upon fix and whose capacity to distinguish between A and quiessence on the 
one hand, and from kinessence on the other, actually increases with redshift (see Fig. 2). The present CMB, SNe and 
galaxy clustering data strongly suggest that a < 1 for quintessence with the inverse power-law potential [7]. However, 
even then the Statefinder remains a useful diagnostic as will be shown below. 

Let us consider another form of kinessence. Below we determine the value of the Statefinder pair for the simplest 
of brane cosmological models - the Dvali-Gabadadze-Porrati (DGP) model [9]. It is important to note that in this 
model 'dark energy' is not the energy associated with a new form of matter, rather its origin is geometrical in nature 
and is entirely due to the fact that general relativity is formulated in 5 dimensional space-time. The model has only 
one adjustable parameter r c - the scale beyond which gravity becomes five-dimensional. This scale can be related to 
the current values of H and VL m0 by the relation H r c = 1/(1 — Q m o)- The FRW equation for this model reads: 



H = J^l + ^ + ±. (11) 

V 3 4r 2 2r c V ' 

(the choice of sign in front of the last term in the right-hand side corresponds to the 'BRANE2' class of models 
according to the terminology of [10]). 

The solution to (11) can be written in the following parametric form: 

. , 2/3 , 3t , 1 - e" 2 ^ 
a = ai sum ' ip , - — = ip -\ , 
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The values of the deceleration parameter and the Statefinder pair read: 

q i + n m ' (i + fi m ) 3 ' (i + ^ m ) 2 ' 1 J 
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In particular, r = 0.74, s = 0.11 for TO = 0.3. At large redshifts the universe becomes matter dominated and r — > 1, 
s -> 0.5. 

At the end of the paper, we estimate the accuracy with which the Statefmder pair (averaged over a range of z) 
can be determined in future SNAP-type satellite missions. The 'SuperNovae Acceleration Probe' (SNAP) is expected 
to observe approximately 2000 type la supernovae within a year up to a rcdshift z <~ 2 and to improve luminosity 
distance statistics by over an order of magnitude [21]. Measurement of the luminosity distance Dl(z) allow us to 
determine the Hubble parameter, since [11,4] 



H(z) 



d (D L {z) 



(14) 



dz \ 1 + z 

To determine the Statefmder pair we use the following model independent parameterization of H(z): 

H 2 (x) = H 2 [a„ x 3 + A + Bx + Cx 2 ] , (15) 

where x = 1 + z and A + B + C = 1 — Cl nl Q. This form is simpler than that used in [12] but it is sufficient for our 
purpose. It becomes exact in the case of the ACDM model (i.e., dark energy being a cosmological constant). Note 
that the fact that we parameterize H 2 (z)/Hq by a 3-parameter fit means that the real H(z) curve is smoothed over 
some redshift interval z ~ z max /3. In principle, the value of f) m0 can be somewhat larger than the current density 
in CDM + baryons if dark energy has a tracker component having equation of state equal to that of matter at high 
z. However, the difference between S7 m0 and f2 m o (if exists at all) is known to be small: f2 m o l-10 TO o. Supernova 
observations of D L and relations (14) and (15) can be used to determine A, B, C and the Statefmder pair {r, s}, since 

(B + Cx)x 

r = 1 — — 



n m x 3 + A + Bx + Cx 2 



2(B + Cx)x 

' 3(3A + 2Bx + Cx 2 )' { ' 

In Fig. 3 we present the results obtained from 1000 random simulations of a SNAP-type experiment for the 'mean 
Statefmder statistic' 



1 f 2 -— 
r = / r(z) dz , 



(17) 



1 f z ™* 

s= / s{z)dz (18) 
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with z max = 1.7. The simulated numbers of SNe la events for one year period of observations are taken to be 50, 
1800, 50 and 15 for the redshift intervals (0 - 0.2), (0.2 - 1.2), (1.2 - 1.4) and (1.4 - 1.7) respectively. The statistical 
uncertainty in the magnitude of SNe is assumed to be constant over redshift and is given by cr mag = 0.15. Details will 
be presented in a companion paper [15]. Fig. 3 shows that a future SNAP-type experiment determining {r, s} can 
easily distinguish a fiducial ACDM model from several alternative time-dependent forms of dark energy, including the 
inverse power-law quintessence potential V cx (f)~ a with a ~ 1 and the DGP brane cosmological model. 
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